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The atomic structure and mechanical properties of the carbyne (monatomic linear chains), con¬ 
taining from 2 to 21 carbon atoms, are theoretically investigated by ab-initio methods. We demon¬ 
strate the existence of a stable cumulene-structure in the inner part of chains with the number of 
atoms N > 10. We present a general stress-strain diagram of chains until the moment when they 
break, which enables to determine their strength, elasticity and fragility. For chains with N ^ 4, 
the relationship between the strength of the chain and the binding energy of the edge atom in 
the chain is established. The existence of scale-effect and ’’even-odd” effect for such properties as 
strength, elasticity and fragility is observed. We demonstrate that the 5-atom carbon chains show 
the maximum strength value. 
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I. INTRODUCTION 

Carbynes (monatomic linear chains of carbon) have 
recently attracted much attention due to their unusual 
physical properties — , and promising applications —AS. 
The potential implementations of these unique functional 
properties essentially depend on the strength and elas¬ 
ticity of the carbine. Moreover, the possibility of obtain¬ 
ing carbine by unraveling out of nanotubes or graphene 
sheets is governed by it mechanical properties Ai£ i n i 12 . 
Therefore, a lot of studies concernig the investigation of 
strength and stability of monatomic carbon chains have 
appeared recently —~ 16 . The results of ab-initio calcu¬ 
lations are still forming the main source of information 
about the mechanical properties of carbyne. Only re¬ 
cently the experimental results of tests to determine the 
tensile strength of carbon atomic chains using a high- 
field method have been published —and extremely 
high level of strength of these chains, which exceeds 270 
GPa was found—. Moreover, their high-field-evaporation 
stability was ascertained. 

At the same time, most of the studies do not present 
the results of targeted investigations of the dependence 
of the strength of chains on their atomic structure. The 
mechanical properties data is usually the auxiliary one 
when analyzing electronic, magnetic and other functional 
properties. Currently, to the best of our knowledge, this 
information is strongly contradictory in literature. For 
example, the data on the strength of monatomic carbon 
chains in various articles differs by the order of magni¬ 
tude—], and in most cases, the properties of chains of 
infinite length are investigated. However, the structure 
and properties of a finite carbon chain (carbyne) differ 
significantly from those for an infinite chain 20 i 21 . The 
present study is focused on ab-initio simulations of car¬ 
byne atomic structures with a detailed analysis of their 
strength, elasticity and fragility. 


II. THE SUBJECT AND THE METHODS OF 
INVESTIGATION 


As it was shown in— iSi, the interatomic distance in the 
chains of different lengths is not constant and depends on 
the number of atoms in the chain. In the present study 
we carry out a detailed analysis of this effect using ab- 
initio calculations of the electronic and atomic structure 
of chains of different lengths with even and odd number 
of atoms. The number of atoms in the chain was varied 
from 2 to 21. 

The mechanical properties and total energies of car¬ 
bon chains were calculated using the quantum-espresso 
program package —, which uses planewaves as a ba¬ 
sis set. We employed ultrasoft pseudopotentials, which 
were custom-generated for carbon atoms using the Van¬ 
derbilt code (v 7.3.4)— The planewave basis cutoff of 
E cu t = 450eU was used in all calculations. To check 
the pseudopotentials and the basis cutoff value, a series 
of test calculations for infinite chains with the structure 
of cumulene and polyyne were performed. The obtained 
interatomic distances and total energies were consistent 
with the results of the previous study—. To calculate and 
visualize the spatial distribution of the electron density in 
the chains, an all-electron Linearised Augmented Plane 
Wave ( LAPW ) method was employed, as implemented in 
the WIEN2K program package— The generalized gradi¬ 
ent approximation (GGA) for exchange-correlation func¬ 
tional was used in all calculations —. 

As the results of ab-initio calculations, the dependence 
of the tensile force F on the value of total chain strain e 
was determined, as well as the values of the bond strain 
e between the first and the second atoms from the edge 
of the chain (Fig. |TJ) . The force acting on the edge atom 
was calculated as: 
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FIG. 1. The dependence of force F on strain e for the whole chain with odd (a) and even (b) number of atoms, and on the 
value of strain e at the edge atomic bond in ’’odd” (c) and ’’even” (d) chains: F c is critical stress for instability of atomic bond 
(bond strength); s c and e c are critical strain of instability of atomic bond and the chain as a whole, respectively; £/ and e/ are 
fracture strain for atomic bond and chain, respectively; fey is coefficient of elasticity of the chain. 


F = 


dE 
da ’ 


(1) 


where E is the total energy of the system; a is the 
current distance between the first and the second atoms 
(length of the edge bond). 

The strain of the whole chain, e, was estimated as: 


e = 1,1 (£) ■ (2) 

where Iq and l are the initial and current chain lengths. 
The strain of interatomic bonds between the first and the 
second atoms was calculated as: 


£ = In 



( 3 ) 


where a o and at are the initial and current length of the 
edge bond, respectively. The force-strain dependencies 
(Fig. [T]) enabled us to estimate the chain strength F c , 
which was calculated as a maximum force value at the 
moment of chain instability. The values of critical strains 
of instability of the edge atom bond with the chain e c 
(Fig. QJ:, d), as well as the value of critical strain for the 
whole chain e c (Fig. [l]a, b) were obtained. The values of 
the elasticity coefficient (stiffness) ky and of the elastic 
Y modulus were calculated as: 


ky 


dF 


de 


e=0 


( 4 ) 



( 5 ) 


where S is the effective cross-sectional area of the 



























3 



FIG. 2. (Color online) The distribution of the electron density 
in the bulk of the chain and its cross-section passing through 
the middle of the edge atomic bond in the unloaded state 
(dashed line) and at the moment of instability of the chain 
(solid line). 

chain, through which the atoms interact with each other. 
To estimate this value, the calculated distribution of elec¬ 
tron density in the chain have been used. Figure [2] shows 
the spatial distribution of the electron density for a chain 
of seven atoms in the unloaded state, calculated by the 
LAPW method. The figure shows the change in the den¬ 
sity of the electron charge in the cross-section of the chain 
in the unloaded state and at the critical strain e c of chain 
instability. The cross-section for the unloaded state was 
chosen in the middle of the distance between the two 
edge atoms. For the critical state of the chain the cross- 
section was chosen at a distance from the edge atom, at 
which the magnitude of the electron charge in a plane 
perpendicular to the axis of the chain, reaches its mini¬ 
mum value. 

This approach enables us to give the lower estimate for 
the transverse size of the area of the force interaction be¬ 
tween the first and the second atoms in the chain, which 
value was found to be d « 2A. We note that the same 
effective diameter of the chain has been used in — to sim¬ 
ulate the formation of carbyne by unraveling. This value 
coincides with the effective magnitude d obtained in 27 
to analyze the thermal conductivity of chains. Further¬ 
more, this value enables to obtain the calculated value 
of strength of carbyne, which correlates reasonably well 
with its experimental value, as established in —. 

III. RESULTS AND DISCUSSION 

Figure [3] shows the distance between the two nearest 
neighboring atoms ai,i+i(i = 1,..., N) in carbon chains 
of different lengths. According to these results, the dis¬ 
tance between the atoms depends on their positions in 
the chain, as well as on the length of the chain itself. This 
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FIG. 3. The dependence of the length of interatomic bond, 
aij, on the number of atoms in the chain N. 

represents the principle difference between the atomic 
structure of carbyne and that of infinite carbon chains. 
The distance between the two outer atoms (the first and 
the second - 011 , 2 ) appeared to be the largest one, while 
the smallest distance was observed between the third and 
the fourth atoms. In chains with less than 16 atoms, 
the distance 01,2 depends on both the total number of 
atoms in the chain, and on whether this number is odd 
or even. Therefore the ’’scale” and ’’even-odd” effects 
occur simultaneously. Our calculations showed that in 
carbynes with more than 10 atoms, the interatomic dis¬ 
tances within the chain (starting from the fourth atom) 
are equal to those in cumulenes, i.e. the internal struc¬ 
ture of carbyne is a cumulene one. It should be noted 
that the cumulene structure is unstable in the chains of 
infinite length, and for such chains a polyyne structure is 
energetically favorable 16 i 28 . Therefore, the existence of 
a cumulene structure in the inner part of the chain is a 
specific characteristic of the atomic structure of carbyne. 
The presence of the edge atoms appears to be the reason 
for stability of a cumulene-structure in the central part 
of the finite length chains. 

Similar to the case of the infinite chains, the ratio 
of the bond lengths in finite chains can be character- 
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FIG. 4. The dependence of the value of bond length alterna¬ 
tion (BLA) on the number of atoms in the chain N: BLAi : 3 
and BLA 2,4 are the values of BLA for atomic bonds a 1,2 and 
02,3; a-2,3 and 03,4 respectively. 


ized using the BLAi j (bond length alternation) quantity. 
The regularities of changes in BLAi^ = 01,2 — 02,3 and 
-BLA. 2,4 = 02,3 — 03,4 values in chains of different lengths 
are presented in Figure [I] General regularity of changes 
in the BLA 13 and BLA. 2,4 values lies in the fact that 
the chains with an even number of atoms, show maxi¬ 
mum values of these quantities, while the chains having 
an odd number of atoms - the minimum ones. Thus, 
the opposite trends in BLA change with the growth of 
the number of atoms N are observed. In even-numbered 
chains the values of BLAi^ and BLA 2 : 4 , decrease with 
N (at N > 6 ), while in odd-numbered chains, they are 
increasing. At N > 19, these values no longer depend on 
the total number / parity of atoms and they approach 
a constant value of BLA ss 0.014A. This value is sig¬ 
nificantly smaller than the BLA for the infinite polyyne 
chain (0.070-0.090A) —. 

Significant changes in the interatomic distances along 
the chain (Fig. [3]) indicate that the interatomic inter¬ 
action energy also varies depending on the position of 
atom in the chain. Therefore, to describe interatomic in¬ 
teraction in the carbyne, the binding energy of each atom 
with the rest of the chain, E b (N), should be estimated. 


The total energy of a finite chain, E(N ), consisting of N 
atoms, can be represented as: 

N 

E(N)=J2E b dN) + NE at , ( 6 ) 

i 

where E at is the energy of a free carbon atom. In car¬ 
byne structures with the number of atoms N > 6 , three 
types of carbon atoms with different binding energies can 
be identified. Our analysis of the interatomic distances 
(Fig. [3]) shows that the first (from the edge) atom should 
have the lowest binding energy, as the distances aq ,2 be¬ 
tween the first and the second atoms are the largest ones 
in carbyne. At the same time, the atoms in the cen¬ 
tral part of the carbyne are expected to have the highest 
values of the binding energies, as the distances between 
those atoms are approximately the same and are close to 
the bond length in cumulenes. The second atoms from 
edge are expected to have intermediate values of the bind¬ 
ing energies. This allows us to represent total energy (6) 
as follows: 

E(N) = NE at + 2 [E b (N) + E b (N)\+(N-4)E b cum , (7) 

where E at is the energy of a free carbon atom, E b (N) 
and E 2 ( N ) are binding energies of first and second atoms 
from chain edge, and E b um = —7.71eF is a binding en¬ 
ergy of the carbon atom in cumulene structure. The total 
energy calculations of finite chains with different num¬ 
bers of atoms allowed us to evaluate the binding energy 
of carbon atoms with different location inside the chain. 
The results of the calculations are presented in Fig. [5] 
In the first approximation, the value of E b ( N ), in chains 
with N > 16 is equal to the average binding energy of 
carbon atom in a chain of three atoms E b = —5.80eV 
(Fig. [5]). Such approximation is based on the fact that 
interatomic distances in the chain are close to the value 
ai 2 in chain with number of atoms N > 16 (Fig. [3]). 
Binding energy of the near-edge atom (E 2 = —6.58eF) 
was obtained from the expression (7) for the total energy 
of the chain containing 16 atoms using the E b = —5.80eF 
value. Interestingly, the average energy of the atom in 
the 5-atom chain is equal to 6.55eV and is close to the 
value of E 2 . This can be explained by the fact that the 
interatomic distances in a 5-atom chain are equal, and 
are close to the value of < 22,3 in the chains with N > 16 
(see Fig. [3]). According to our calculations, the value 
of E b , depends both on the total number of atoms in 
carbyne and on the parity of this number (see Fig. [5]). 
In chains with an odd number of atoms, the values of 
E b are greater in magnitude and decrease with increas¬ 
ing number of atoms in the chain, approaching the value 
of E b = —5.80eF. The situation is different for even- 
numbered chains, where binding energy of the edge atoms 
increase with N (see Fig. [5]). It is interesting to note the 
fact that the stronger odd-numbered carbynes are insu¬ 
lators, and even-numbered ones are conducting systems. 
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FIG. 5. The dependence of the binding energy on the number 
of atoms: E b (N) and E%(N) - for the first and second atoms 
from chain edge, E b um - for inner atoms, A 6 (TV) - energy 
of separation, calculated on the equation (8), Eq(N) - the 
average binding energy per atom (solid line - calculated on 
the equation (11)). 


However, a detailed study of the electronic structure of 
carbyne systems is beyond the scope of the present paper 
and should be a topic of a separate investigation. 

The data presented in Fig. [5] demonstrates that the 
atomic bonds of the edge atoms are the weakest ones. As 
we show below, this fact leads to the edge atoms tearing 
off upon the stress application. Formally, the energy of 
separation of one atom from a chain can be written as 
follows: 


E b (N ) = E(N) - E(N - 1) - E aU (8) 

where E(N) and E(N — 1) are the total energies of 
chains with TV and TV— 1 atoms respectively. The regular¬ 
ities of change in this value with increasing of the number 
of atoms in carbyne are presented in Fig. [5] By substitut¬ 
ing E(N) and E(N — 1) from (7) to (8), while taking into 
account that E^N) = const , one can see that the non¬ 
monotonic change in E b (N ) value is mainly due to the 
change in the binding energy of the edge atom, E b (N ): 


E b (N) = 2 [E b (N) - E b (N - 1)] + E b cum , (9) 

At TV > 16 the E b (N) value is almost equal to 
E b (N — 1). Therefore, with increasing of TV, the value of 
E b (N) approaches the value of binding energy of the car¬ 
bon atom in cumulene, E cum = —7.71eV, which agrees 
well with the results of ab-initio calculations (Fig. [5]). 
Therefore, the equation (8) describes the energy balance 
in the separation of one atom from the chain, so it can 
be used to determine the binding energy in the infinite 
chain. However, it does not take into account specific 
features of interatomic interaction in the chains of finite 
length. 

Equation (7) for the total energy of chain allows us to 
write down analytical expression for the average binding 
energy Eg (TV) per atom: 


E b 0 (N) 


2[E b (N)+E b (N) 

TV 


2 Et 


+ E 


b 

cumi 


( 10 ) 


When TV > 16 the values E b (N) and E^N) do not 
depend on the number of atoms in a chain, and are equal 
to E b (N) = -5.80eF and E b (N) = -6.58eH. In this 
case the average value of binding energy per atom can be 
found using the following expression: 


E b (N) = ^+E 


b 

cum 5 


(ii) 


where A = 6.08eH. 

According to Fig. [5] this relation agrees well with the 
results of our ab-initio calculations and can be used to 
predict the average value of the binding energy per atom 
in carbynes containing different number of atoms. The 
difference between E(N) for the chain of finite length and 
infinite chain with a cumulative-structure is due to the 
presence of the edges atoms (edge effect), and in case of 
N = oo we have E b (N) = E b cum . 

To calculate the mechanical properties of carbyne, 
we performed modeling of tension of chains of various 
lengths up to their complete break. As results of these 
calculations, the values of elasticity coefficient fcy and the 
elasticity modulus Y, the value of the maximum force F c 
as well as the respective values of critical strains of the 
whole chain e c (Fig. [ljj, b) and the edge bond s c (Fig. |T] 
c, d) were determined. Besides these values, the breaking 
strain of the whole chain e/ (Fig [I] a, b) and the edge 
bond e/(Fig. |T] c, d) were found. The F c value defines 
the level of strength of carbine, and e c characterizes its 
fragility. 

According to the results of our ab-initio simulations, 
the magnitude of strength of the carbyne is determined 
by the strength of atomic bond of the edge atom. The 
value of this strength depends both on the total num¬ 
ber of atoms in the chain, and on whether this number 
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FIG. 6. The effect of the number of atoms in carbyne on its 
strength F c , fragility e c , and hardness ky- 


Binding energy, E (eV) 

FIG. 7. The relationship between the strength of the chain, 
F c and the binding energy Ei of the edge atom. 


is odd or even. We observe that the strength of car- 
bynes with an odd number of atoms is higher in com¬ 
parison with the strength of even-numbered carbynes. 
Carbyne containing 5 atoms has maximum strength of 
F c = 13.09nN. Chains with an odd number of atoms 
not only show higher strength, but are also less fragile 
because their instability takes place at higher values of 
critical strain e c . The simulation results, presented in 
Fig. [6l show that the difference between the strength of 
even-numbered and odd-numbered chains decreases with 
increase of the number of atoms in chain, and at N > 12 
this difference practically disappears. The increase in the 
number of atoms in the chains gives rise to an increase 
in the stiffness of the chain, ky , and to a steep (1.3 - 1.5 
times) increase in their fragility (decrease in e c ). The lat¬ 
ter is due to the inhomogeneous deformation of the chain, 
which shows itself as the localization of deformation in 
the bond between the first and second atoms because the 
edge bond is weakest. This effect is enhanced with in¬ 
creasing the chain length. Quantitatively, it shows itself 
in a significant difference between the value of fracture 
strain of the whole chain , e/, and the critical fracture 
strain of the atomic bond between the first and second 
atoms Ef (Fig. [TJ. 

It should be noted that the disappearance of the 
’’scale” and the ’’even-odd” effects for strength occurs at 
reaching of a smaller number of atoms in the chain (12 
atoms) instead of 16 atoms for the interatomic interaction 
energy. For the coefficient of elasticity, ky, the difference 


between ’’even” and ’’odd” chains vanishes starting from 
N > 16 atoms. 

Usually the strength of interatomic bonds is estimated 
not directly, but using the value of binding energy. 
Therefore, we did a comparison between the binding en¬ 
ergy of the edge atom E\, and the value of the critical 
force of instability of the edge bond, F c , which defines the 
strength of the whole chain. As it is shown in Fig. [7] 
for chains containing more than 3 atoms, the increase in 
the binding energy of the edge atom, E\, is accompanied 
by the increase in strength of the chain. 

According to our results, the strength of carbyne 
ranges from 11.3 nN (360 GPa) to 13.1 nN (417 GPa). 
These values are higher than the experimental strength 
of 270 GPa—. It is probably due to the fact that the ex¬ 
perimental value of 270 GPa is a lower bound of strength 
of carbyne, as it determines the strength of the contact 
atomic bond between the chain and graphene from which 
this chain was drawn out. 

Finally we should note that the different symmetry 
of the even-numbered and odd-numbered carbyne is the 
main reason for ’’even-odd” effect. One carbon atom is 
always located in the center of symmetry of the chain 
with an odd number of atoms. This results in the fact 
that the binding energy of such atom is the highest one. 
By virtue of the long-range interatomic interaction in the 
carbon chains, it causes an increase in bond strength of 
the edge atoms. 
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IV. CONCLUSION 

We performed a detailed theoretical study of atomic 
and mechanical properties of the carbyne chains using 
the ab initio methods. Our results demonstrate that in 
carbynes with number of atoms TV > 10 the length of the 
” inner” interatomic bonds (starting from the fourth from 
the edge) becomes equal to the bond length in cumulenes. 
The existence of this inner structure is one of the most 
fundamental differences between the atomic structure of 
carbyne and that of the infinite chains, which show stable 
polyyne structure. 

It is shown that unlike the carbon chains of infinite 
length, the interatomic distance in carbynes varies along 
the chain : The bond length between the first and the 
second (from the edge) atoms is the largest one, while the 
bond length between the third and the forth atoms has 
the smallest value. These atoms are effectively separating 
the inner part of carbyne from the first and second atoms. 
In short carbines (TV < 16), the distance between the first 
and second atoms depends both on the total number (TV) 
of atoms in the carbyne, and on whether N is odd or even. 
With the growing number of atoms in carbyne, the length 
of the edge interatomic bond in even-numbered carbynes 
decreases, while this bond is increasing in odd-numbered 
carbines. When the total number of carbon atoms is 
TV > 16 , these values approach 1.306A, which is only 
0.005A greater than the bond length in polyyne. 

The total energy calculation showed that the energy of 
interatomic interaction in carbynes depends on the posi¬ 


tion of the atom in the chain. Depending on the value 
of the binding energy, one can distinguish three types of 
atoms: (i) the edge atoms that have the lowest binding 
energy (=-5.80eV for TV > 16); (ii) inner atoms which 
binding energy in chains with TV > 10 is equal to the 
binding energy in cumulenes (-7.71 eV); and (iii) the in¬ 
termediate atoms with a binding energy of -6.58 eV. In 
short carbynes (TV < 16) the binding energy of the edge 
atoms depends on both the total number of atoms in the 
carbyne, and the parity of this number. The binding en¬ 
ergy of these atoms in the even-numbered carbynes is al¬ 
ways lower. Its value increases with the number of atoms, 
and in the odd-numbered carbynes, it, on the contrary, 
decreases. 

We show that the mechanical properties of carbyne 
containing more than 4 and less than 12 atoms are gov¬ 
erned by the scale and ’’even-odd” effects. The chains 
with an odd number of atoms are stronger and less fragile 
than the even-numbered chains. However, these effects 
vanish when the number of atoms in the chain is more 
than 10. 

Carbyne is the strongest material in the world. Lower 
experimental estimation of its strength is equal to 270 
GPa at 3K. Our ab initio calculations demonstrate that 
the 5-atom carbine chains show the maximum strength 
of 13.1 nN (417 GPa), while 392 GPa was obtained for 
this value for the the carbynes with TV > 12. 
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